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ABSTRACT 

A single  source  transportation  problem  is  an  ordinary  transportation 
problem  with  the  additional  requirement  that  the  entire  demand  at  each 
demand  location  be  supplied  from  a single  supply  location.  It  is  a special 
case  of  Ross  and  Soland's  generalized  assignment  problem.  Such  problems 
occur  frequently  in  applications.  This  paper  gives  two  heuristic  solution 
methods  and  a branch  and  bound  algorithm  for  solving  single  source  trans- 
portation problems.  A discussion  of  the  branching  rules,  variable  fixing 
rules,  and  the  computation  of  weak  lower  bounds  is  given.  Computational 
experience  with  the  solution  of  randomly  generated  problems  having  up  to 
40,000  integer  variables  is  reported. 


1 . INTRODUCTION 

In  this  paper  we  consider  ordinary  transportation  problems  with  the 
additional  restriction  that  each  demand  must  be  entirely  supplied  from  a 
single  source.  It  is  therefore  a special  kind  of  generalized  assignment 
problem  in  the  sense  of  Ross  and  Soland  [ 5 ] . 

There  are  many  applications  in  which  such  requirements  are  made  on 
the  solution.  For  instance,  the  supplying  of  supermarket  orders  from  a 
network  of  central  warehouses  frequently  has  this  restriction.  In  military 
applications  it  is  common  to  require  that  all  troops  going  on  the  same 
mission  leave  from  the  same  staging  area.  When  a group  of  computers  is  used 
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to  fill  a set  of  computation  demands,  the  common  requirement  is  that  all 
computation  on  each  single  job  be  performed  entirely  by  a single  computer  [1]. 
Many  facility  location  models  contain  a single  source  requirement  on  the 
shipment  of  goods  from  the  opened  facilities  to  the  demand  locations  [6] . 

Single  source  transportation  problems  are  also  related  to  a large  class  of 
"loading"  or  "packing"  type  problems  where  we  attempt  to  assign  a set  of 
weighted  objects  to  boxes  or  bins  which  have  weight  capacities  [2  ] . 

In  [ 3 ] De  Maio  and  Roveda  formulated  a problem  which  is  a special 
version  of  the  generalized  assignment  problem  stated  later  by  Ross  and 
Soland  [5],  Srinivasan  and  Thompson  [ 7]  showed  how  to  transform  De  Maio 
and  Roveda’ s problem  into  a single  source  transportation  problem.  They 
proposed  solving  the  latter  using  a branch  and  bound  cost  operator  algorithm 
which  used  the  ordinary  transportation  problem  as  a relaxation  of  the  single 
source  problem. 

In  the  present  paper  we  describe  and  give  computational  results  for  two 
heuristic  solution  methods  and  a cost  operator  algorithm  which  is  similar  to 
the  algorithm  described  in  [ 7 ] . The  present  algorithm  differs  from  that  in 
[7]  in  the  following  respects:  (a)  We  have  replaced  the  "row  unique" 

solution  concept  in  [ 7 ] by  our  "single  source"  solution  concept,  see  section  5; 
(b)  Different  variable  selection  and  branching  rules  are  used.  (c)  Weak 
lower  bounds  are  calculated  and  are  used  in  fathoming  as  well  as  for  variable 
selection.  (d)  A non-basic  variable  fixing  rule  has  been  added. 

In  section  6 computational  results  from  the  solution  of  problems  ranging 
in  size  from  5 x 10  to  100  by  400  are  presented.  All  of  these  problems  were 
generated  randomly  using  the  method  described  in  Ross  and  Soland  [ 5 ] . A 
discussion  is  given  concerning  the  way  that  problem  difficulty  depends  on  the 
setting  of  parameters  in  the  Ross  and  Soland  problem  generator. 
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2.  STATEMENT  OF  THE  PROBLEM 

In  the  single  source  transportation  problem  we  consider  a set  of 
sources  I = each  having  capacities  a^  > 0;  a set  of  uses  (or 

users)  J = {l,...,n}  each  having  known  demands  > 0;  and  a set  of 
costs  c^  of  supplying  use  j from  source  i.  The  problem  is  to  assign 
sources  to  uses  so  that:  (i)  the  total  amount  shipped  from  each  source  does 
not  exceed  its  capacity;  (ii)  each  use  is  supplied  by  exactly  one  source; 
and  (iii)  the  total  cost  Z of  the  assignment  is  minimized. 

We  shall  assume,  as  a necessary  but  not  sufficient  condition,  for 
feasibility  that 

S a > S b (1) 

iel  1 jej  J 

i.e.,  that  supply  exceeds  or  equals  demand.  Defining  j'  » J U (n+lj , 

c.  . » 0 for  iel,  b , « £ a.  - £ b.,  and  letting  x. . be  the 

i,n+l  n+1  i8].  i jeJ  j ij 

amount  shipped  from  i to  j,  we  can  write  the  single  source  transportation 
problem  as: 


Minimize  Z - £ £ c x 

iel  jeJ  J J 


Subject  to 


£ ,x.  • a for  iel 

jeJ  J 1 


£ x.  - b for  jej' 


ij  J 


x^j  > 0 for  iel,  jej 
x^  * 0 or  bj  for  jeJ 
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The  problem  defined  by  (2)- (5)  is  an  ordinary  transportation  problem 
which  we  will  call  problem  P.  We  denote  by  p'  the  problem  in  (2)- (6). 

We  will  say  that  problem  P is  the  transportation  relaxation  of  problem  p'. 

REMARK:  In  the  case  where  all  of  the  demands,  bj  jcj,  are  equal, 

?'  can  be  transformed  into  an  equivalent  (semi-assignment)  problem,  P*, 
having  integral  supplies  and  all  demands  equal  to  one.  This  can  be  done  by 


making  the  following  transformations 

i:  let  bj  ■ b for  all  JcJ 

★ 

Cij  " Cij  X b 

id  j«J; 

(7) 

* 

ai  * di 

iel 

(8) 

where  ■ d^  X b + k^ 

0 < kt  < b 

(9) 

★ 

bJ  ” 1 

JeJ 

(10) 

* * 

bn  ,i  - E a - n 
n+1  id  1 

• 

(11) 

The  reader  can  verify  that  by  making  the  transformations  (7)  - (11)  the 

★ 

resulting  single  source  transportation  problem,  P , is  equivalent  to  the 

/ * * 
original  problem,  P , in  the  sense  that  X is  a solution  to  P if  and 

only  if  X7  - bX  is  a solution  to  p'.  The  significance  of  transforming  p' 

* * 

into  P is  that  P is  a transportation  problem  with  integral  supplies  and 
unit  demands.  Such  problems  are  known  as  semi-assignment  problems.  Thus  due 
to  the  well  known  unimodularity  property  of  the  basis  matrix  for  a transporta- 
tion problem,  we  know  that  the  solutions  to  the  transportation  relaxation  of  P 
will  either  be  0 or  1 and  therefore  will  satisfy  the  single  source  require- 
ment automatically  without  any  special  search  algorithm.  In  this  paper  we  con- 
sider the  case  of  unequal  demands. 
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DEFINITION  1.  By  a single  source  basic  solution  to  the  transportation 
problem  P we  mean  a feasible  basic  solution  with  the  property  that  for  each 
j«J  there  is  a row  index  i ■ i(j)  such  that  x^  ■ b^;  in  other  words  a 
solution  in  which  each  demand  is  completely  supplied  by  a single  source. 

In  [ 7 ] the  idea  of  row  unique  solutions  were  introduced,  which  are  a 
special  kind  of  single-source  solutions.  A row  unique  basic  solution  to 
problem  P is  a basic  feasible  solution  with  the  property  that  for  each  jeJ 
there  exists  a unique  row  i ' - i(j)  such  that  xj/j  is  the  only  basic 
variable  in  column  j.  For  a nondegenerate  transportation  problem,  these  two 
concepts  are  identical.  In  the  case  of  (primal)  degenerate  problems,  it  is 
necessary  to  use  the  single-source  solution  concept  instead  of  the  row  unique 
concept.  We  will  elaborate  on  this  point  in  Section  5.  Since  many  of  our 
problems  are  degenerate  we  concentrate  on  the  former  concept  here. 

Any  single  source  solution  to  the  transportation  problem  P gives  rise 
to  a feasible  solution  to  problem  p'.  In  Section  4 we  describe  a branch  and 
bound  algorithm  which  solves  p'  by  finding  single-source  solutions  to  a series 
of  transportation  problem  relaxations  of  p'.  Each  of  these  relaxed  problems 
differs  from  P by  having  the  cost  of  some  of  the  cells  set  equal  to  +M  or  -M 
(where  M stands  for  a number  much  larger  than  the  absolute  value  of  any  cost) 
in  order  that  those  cells  are  forced  into  or  out  of  the  basic  solution  for  the 
corresponding  relaxed  problem. 

3.  HEURISTIC  SOLUTION  METHODS 


In  order  to  reduce  the  size  of  the  search  tree  in  the  branch  and 
bound  search  process  we  developed  two  heuristic  solution  methods  that  almost 
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invariably  find  feasible  solutions  to  p7  in  a short  time.  The  smallest 
value  of  the  objective  function  for  such  heuristic  solutions  is  used  as  an 
initial  upper  bound  in  the  branch  and  bound  algorithm.  In  many  cases  these 
heuristic  methods  actually  find  an  optimal  solution,  as  will  be  discussed  in 
Section  6 where  data  from  problem  solutions  is  given. 

The  two  heuristic  solution  methods  we  used  differ  only  in  the  order 
in  which  uses  (columns)  were  selected  to  be  assigned  to  sources  (rows).  The 
first  method  calculates  "regrets"  similar  to  those  of  the  VAM  starting 
solution  method  [4]  for  transportation  problems;  the  second  method  selects 
uses  in  the  order  of  non-increasing  demand  sizes. 

To  describe  the  regret  heuristic  let  c . be  the  kth  smallest 

JkJ 

cost  in  column  j.  Define  Reg(j),  the  regret  for  use  j to  be 


Reg(j)  - <c  , - c )b  (12) 

J2J  J 

The  first  use  to  be  assigned  to  a source  is  one  whose  regret  is  largest.  Once 
it  has  been  assigned  (in  a manner  to  be  discussed  in  the  next  paragraph),  it 
is  removed  from  the  set  U of  unassigned  uses,  new  regrets  are  calculated 
as  necessary,  and  a second  use  with  largest  regret  is  chosen  to  be  assigned,  etc. 

Given  that  use  j is  to  be  assigned,  we  first  construct  a set 
of  sources  to  which  it  can  be  assigned,  according  to  the  following  rule: 


SJ  ‘ tk,l|ckJ  ' S'j1 


where  is  a parameter  chosen  by  some  rule  such  as 


(13) 


a 


J 


(14) 


In  other  words  S.  consists  of  the  indices  of  all  sources  whose  costs  in 
J 

column  j differ  from  the  smallest  cost  in  column  j by  an  amount  less  than 
or  equal  to  a. . Then  one  of  the  indices  ieS^  is  chosen  randomly  to  be 
the  actual  source  to  supply  demand  b^ . 

The  reasons  for  the  random  choice  among  indices  in  the  set  are: 

first,  with  a specific  choice  rule  it  is  quite  easy  to  make  early  choices  which 
lead  to  infeasible  solutions;  and  second,  with  the  random  choice  rule  we  can 
repeat  the  heuristic  choice  rule  several  times  and  retain  the  smallest  cost 
solution  found. 

In  order  to  state  the  regret  heuristic  more  precisely  we  first  define 
the  notation  to  be  used. 

U = set  of  unassigned  uses 

Z ■ the  cost  of  current  solution  (or  partial  solution) 

S.  ■ set  of  all  sources  kel  such  that  c.  . - c.  . < a.. 

j kj  jjj  - j 

0^  ■ a parameter  whose  size  determines  the  number  of  elements 

in  Sj.  (In  some  problems  it  may  be  desirable  to  have 
Oj  change  as  more  steps  are  taken  in  the  algorithm.) 


(HI)  Regret  heuristic. 

(1)  (Initialization)  Let  U ■ J,  Z “ 0.  For  each  jeU  find 

h’  j2’  and  V Set  CTj  " cj3j  ■ cjlj* 

(2)  (Check  feasibility)  If  c.  , ■ M for  some  jeU  go  to  (8). 

J1J 

Else  go  to  (3). 

(3)  (Choose  the  use  having  the  largest  regret).  Choose  jeU  such 

that  Reg(j)  » (c  . - c ,)b  is  a maximum. 

J2-J  J 
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(4)  (Determine  the  choice  set.)  Calculate  Sj  which  is  defined  in  (12). 

(5)  (Select  a source.)  Choose  ieS^.  at  random.  Make  the  following 
replacements : 

U by  U - {j} 

Z by  Z + Cljb. 

a£  by  - bj  . 

(6)  If  U » <t  go  to  (9).  Else  go  to  (7). 

(7)  (Update  the  costs.)  For  all  heU  if  b.  > a.  set  c.,  ■ M. 

n i ih 

Find  the  two  lowest  cost  cells  c,  , , c.  in  each  column. 

h^n  h^n 

Go  to  2. 

(8)  Current  solution  is  infeasible.  Stop. 

(9)  Feasible  solution  found.  Stop. 

The  second  heuristic,  which  we  call  the  "largest  demand  heuristic," 
is  identical  to  the  regret  heuristic  except  that  in  step  (3),  instead  of 
choosing  at  each  step  the  use  with  the  next  largest  regret,  we  choose  the  use  with 


the  next  largest  demand.  Thus  in  the  largest  demand  heuristic  (H2),  step  3 
is  replaced  by: 

(3 1 ) (Choose  the  next  use.)  Let  jeU  be  the  index  of  an  unassigned 
use  whose  demand,  b , is  largest. 

• J 

Since  step  (3)  requires  more  effort  in  the  regret  heuristic  than  in  the 
largest  demand  heuristic  we  find,  as  expected,  that  the  regret  heuristic  uses 
more  CPU  time.  However,  neither  of  the  heuristics  requires  very  much  time. 

For  example,  five  runs  of  the  regret  heuristic  on  a 75  x 200  problem  requires 
about  4 seconds  of  CPU  time  (DEC  20  computer). 

The  regret  heuristic,  being  the  "greedier"  heuristic,  tends  to  generate 
lower  cost  feasible  solutions  than  the  largest  demand  heuristic.  The  largest 
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demand  heuristic  was  designed  to  find  good  feasible  solutions  in  the  case 
when  the  regret  heuristic  failed  to  find  any  solution.  In  all  of  our  compu- 
tational experience  this  phenomenon  occurred  only  once;  yet  it  may  be  more 
frequent  in  problems  with  a different  data  structure.  We  might  add  at  this 
point  that  neither  of  the  heuristics  can  be  guaranteed  to  find  a feasible 
solution  to  a given  problem.  It  can  be  easily  verified  that  the  problem  of 
finding  a feasible  solution  to  p'  is  NP  complete.  Finding  a feasible 
solution  to  ?'  is  equivalent  to  determining  whether  one  can  find  a partition, 

S, , of  a set  of  integers  {b,,...,b  } (the  demands)  such  that  the 
i m in 

sum  of  the  elements  in  equals  ai  (the  supply)  for  i = l,...,m. 

The  latter  problem  is  known  to  be  NP  complete.  Thus  it  is  probably  necessary 
to  carry  out  a partial  enumeration  of  assignments  such  as  is  done  in  the 
method  of  section  4,  in  order  to  guarantee  the  finding  of  a feasible  solution 
to  P*,  when  such  a solution  in  fact  exists. 

4.  THE  COST  OPERATOR  ALGORITHM 

We  now  describe  a branch  and  bound  algorithm  which  uses  cost  operators 
in  the  sense  of  Srinivasan  and  Thompson  [ 8 ] . The  basic  idea  in  the  algorithm 
is  that  the  solution  of  problem  P*  is  obtained  by  solving  a sequence  of 
related  transportation  problems  P,  Pj_,...,Pt  until  we  have  found  a problem, 
say  P^,  such  that  the  optimal  transportation  solution  to  P^  is  also  an 
optimal  single  source  solution  to  . Each  successive  transportation  problem 
differs  from  its  predecessor  in  that  certain  costs  have  been  changed  from 
their  original  values  to  either  -M  or  +M.  When  we  drive  a cost  to  -M 
we  will  say  we  have  "fixed  in"  the  corresponding  cell,  and  when  we  drive  it  to 
+M  we  will  say  that  we  have  "fixed  out"  the  cell.  We  now  give  more  precise 


definitions . 
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We  shall  say  that  cell  (i,j)  has  been  fixed  out  of  vhe  basis  when  a 
cost  operator  has  been  applied  to  the  problem  so  that  c_  becomes  equal  to  -HM, 
where  M is  so  large  that  x^  = 0 for  all  optimal  solutions  to  the  new  problem. 

The  operation  of  freeing  a cell  (i,j)  which  has  been  fixed  out  is  the  applica- 
tion of  a cell  operator  to  drive  c^  back  to  its  original  value. 

We  shall  say  that  cell  (i,j)  has  been  fixed  in  the  basis  when 

(a)  a cost  operator  has  been  applied  to  the  problem  so  that 

the  cost  c„  becomes  equal  to  -M,  where  M is  so  large 
that  x^  = b ^ in  any  optimal  solution  to  the  new  sm; 

(b)  a^  has  been  replaced  by  a^^  - b^; 

(c)  for  any  k such  that  b^  > a^  we  fix  cell  (i,k)  out  of  the 
basis,  since  source  i cannot  supply  the  demand  at  k. 

By  freeing  a cell  (i,j)  which  has  been  fixed  in  we  mean  to  undo  the  actions 
listed  above,  so  that  the  costs  and  rims  go  back  to  their  previous  values. 

In  (A)  - (D)  below  we  describe  the  various  steps  of  the  algorithm. 

The  algorithm  is  stated  in  detail  in  (E)  and  a simple  example  is  worked  in  (F). 
(A)  Tree  Search  Rules 

The  branch  and  bound  algorithm  which  we  are  going  to  describe  is  of  the 
LIFO  (last  in  first  out)  or  depth  first  variety.  The  method  starts  by  first 
solving  the  problem  P.  If  its  solution  satisfies  the  single  source  property 
we  are  finished.  If  not,  we  apply  a nonbasic  variable  fixing  rule,  to  be 
described  later,  which  fixes  some  of  the  nonbasic  variables  at  zero.  Next 
a cell  is  selected  by  one  of  the  cell  selection  rules  given  later,  and  the  cell 
is  fixed  in  the  basis.  This  procedure  is  continued  until  either  (a)  a single 
source  solution  is  obtained  or  (b)  the  objective  function  value  of  the  new 
problem  exceeds  the  current  upper  bound.  In  case  (a)  we  compare  the  solution 


-11- 


with  the  current  best  single  source  solution  and  update  the  latter  and  the 
current  upper  bound  if  the  current  solution  is  better.  In  either  case  we  back- 
track by  freeing  the  last  cell  fixed  in.  Then  (using  the  LIFO  rule)  we  choose 
the  last  cell  fixed  in  and  consider  it  for  fixing  out.  If  the  objective 
function  value  of  the  current  relaxation  plus  a weak  lower  bound  (to  be 
described  later)  exceeds  the  current  upper  bound,  then  we  continue  backtrack- 
ing. Otherwise  the  cell  is  fixed  out  and  the  search  process  continues  deeper 
in  the  search  tree.  A more  precise  description  will  be  given  later  in  this 
section. 

(B)  Nonbasic  Variable  Fixing  Rule 

At  the  initial  node  of  the  search  tree  the  transportation  problem  P 

is  solved.  If  the  solution  to  P is  not  a single  source  solution  then  it 

is  possible  to  examine  the  nonbasic  variables  in  the  solution  to  P in  order 

to  fix  some  of  them  at  zero.  Let  Z , Z be  the  objective  function  value  of  P, 

p u 

and  the  current  upper  bound  respectively.  Let  u^  and  v^  be  optimal  dual 
variables  associated  with  the  current  basic  solution  to  P.  For  any  nonbasic 
variable  (i,j),  iel,  jeJ,  if 


(c. . -(u.  +v.))*b.  >Z  -Z  -1 
ij  i J j ~ u P 


then  we  can  set  x_  = 0 for  the  remainder  of  the  procedure.  This  is  true 

since  if  a nonbasic  variable  x„  were  to  become  positive  then  it  can  only 

assume  the  value  b^ , in  which  case  the  left  hand  side  of  (15)  represents 

the  minimum  amount  by  which  the  objective  function  value  of  P will  increase 

if  we  were  to  require  that  x. . = b..  Thus  if  this  amount  exceeds  the  current 

ij  J 

gap  between  the  upper  bound  and  the  transportation  relaxation  then  we  can  be 
sure  that  x.^  will  not  be  positive  in  any  single  source  solution  to  ?' . 
Actually  this  variable  fixing  test  could  be  repeated  at  any  node  of  the  search 
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Cree,  in  which  case  Z^  would  be  replaced  by  the  objective  function  value  of 
the  current  relaxation.  However,  the  computational  cost  of  doing  so  would 
probably  exceed  the  benefit. 

It  is  interesting  to  note  that  the  left  hand  side  of  (15)  increases 
as  the  demand,  b^,  increases.  Intuitively  this  seems  plausible  since  forcing 
a larger  amount  of  demand  into  a single  cell  causes  more  of  a restriction  on 
the  problem.  In  other  words  the  constraint  x„  > b^  becomes  stronger  as  b^ 
increases.  We  have  found  that  the  variable  fixing  rule  (15)  will  set  to  zero 
as  many  as  957.  of  the  nonbasic  variables.  When  solving  problems  in  which  60% 
or  more  of  the  nonbasic  variables  are  fixed  to  0 by  the  variable  fixing  rule, 
it  would  probably  be  worthwhile  to  use  a code  which  takes  advantage  of  the 
sparse  structure  of  the  resulting  transportation  problem.  We  have  not  so  far 
used  such  a sparse  code  in  our  studies. 

(C)  Weak  Lower  Bounds 

The  algorithm  that  we  have  implemented  is  essentially  a linear  pro- 
gramming based  implicit  enumeration  scheme.  At  any  node  of  the  search  tree 
we  have  solved  a relaxation  of  problem  ?'  where  a variable,  x^^  , is 
either  fixed  at  zero,  fixed  at  b ^ , or  free  to  assume  any  value  between 
zero  and  b^.  If  the  objective  function  value  of  the  current  problem,  Z c> 
does  not  exceed  the  current  upper  bound,  Z , and  its  solution  is  not  single 
source,  then  we  choose  one  of  the  free  variables,  say  x^ , and  we  create 
two  new  problems.  In  one  of  the  new  problems  we  require  that  x^  ■ 0 and  in 
the  other  problem  we  require  that  x_  =*  b ^ . let  us  denote  their  objective 
function  values  by  Z.  and  Z,  respectively.  By  calculating  weak  lower 


bounds  we  can  determine  the  minimum  amounts  by  which  Z must  increase  if  we 

J c 

require  that  x..  » 0 or  x. . = b,.  Thus  for  any  free  variable  x. we  can 
ij  ij  J ij 

calculate  a weak  lower  bound  WLB_  in  the  x, . * 0 direction  and  WLB,  in 

0 ij  b 

the  x^  ■ b^  direction  where, 
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or 


Z < Z + WLB-  < Zn 
c — c 0-0 


Z < Z + WLB,  < Z, 
c — c b — b 


What  we  hope  Co  find  is  Chat  either 


(16) 


(17) 


Z + WLB-  > Z (18) 

c 0 — u 

or 

Z + WLB.  > Z (19) 

c b — u 

or  both.  In  case  (18)  ((19))  we  can  fathom  the  search  tree  in  a particular 
direction  without  having  to  calculate  Zg  (Z^ > . Of  course  there  is  a cost- 
benefit  tradeoff  involved  here  since  some  computational  effort  is  required 
to  calculate  the  weak  lower  bounds.  Thus  the  real  question  is  whether  the 
extra  information  obtained  from  calculating  the  weak  lower  bounds  is  worth 
the  added  computational  effort.  We  have  found  that  in  almost  all  of  the 
tests  which  we  have  performed  the  weak  lower  bounds  have  proven  to  be  cost 
effective.  The  weak  lower  bounds  are  useful  not  only  for  fathoming  nodes 
in  the  search  tree,  but  for  choosing  a "good"  variable  for  branching.  We 
will  describe  the  branching  rules  which  we  have  selected  in  the  next  section. 


We  now  describe  how  the  weak  lower  bounds  are  calculated.  Suppose  that 
at  a given  search  tree  node  we  wish  to  calculate  a weak  lower  bound,  WLBg, 
in  the  direction  x^  ■ 0.  This  is  done  by  applyinp  a positive  cell  cost 
operator  6c^\  t0  t*ie  ceH  (i,j)  which  determines  the  maximum  amount,  |i+,  by  which 
the  cost  in  cell  (i,j)  can  be  increased  while  maintaining  optimality  in  the  cur- 
rent basis.  This  is  called  the  positive  basis  preserving  cell  cost  operator,  and 
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is  described  in  [ 8 ] . At  this  point  cell  (i,j)  contains  some  flow,  x^  , 
between  zero  and  b^  and  the  weak  lower  bound  WLBQ  is, 

WLBq  - H+  x_  (20) 

WLBg  is  a lower  bound  on  amount  by  which  the  objective  function  value  of  the 
current  relaxation  will  increase  if  we  branch  in  the  direction  x^  * 0. 

For  a proof  see  [ 9 ] • 

In  a similar  but  slightly  more  difficult  manner  we  can  calculate  a weak 
lower  bound,  WLB^*  in  the  direction  where  x„  » b^.  In  this  case  we  apply  a 
negative  cost  cell  operator  6c^j  t0  cell  (i,j)  which  determines  the  maximum 
amount,  p>  > by  which  the  cost  in  cell  (i,j)  can  be  decreased  before  the 
current  basis  becomes  non-optimal.  If  the  cost  in  cell  (i,j)  is  changed  to 
c^j  where 

CIj  " Cij  ' (21> 

then  some  nonbasic  cell,  say  (p,q),  will  have  a reduced  cost,  c - (u  + v ), 

pq  p q 

of  zero.  If  cell  (p,q)  is  added  to  the  basis  then  a unique  "cycle"  or  "loop"  is 
formed  in  the  current  basis  tree.  From  this  cycle  we  can  determine  a cell  (e,f) 
which  would  leave  the  basis,  and  the  amount  A,  by  which  the  flow  along  this 
cycle  would  change  if  cell  (p,q)  were  brought  into  and  cell  (e,f)  were  removed 
from  the  basis.  The  valid  weak  lower  bound  WLB,  is  (see  [9]) 

D 


WLB,  - p.  A 

D 


Notice  that  WLB^  is  more  difficult  to  obtain  than  WLBq  since  in  the  former 
case  we  must  find  the  minimum  giver  amount,  A,  along  the  cycle  created  by 
an  incoming  nonbasic  cell. 


w 


» 
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The  effectiveness  of  the  weak  lower  bounds  depends  for  the  most 
part  on  how  close  the  current  upper  bound,  Z^,  is  to  the  value  of  the 
current  relaxation,  Z^.  If  Z^  - Z is  small  then  (18)  and  (19)  provide 
strong  fathoming  devices.  Thus  the  ability  of  the  heuristics  to  generate 
low  cost  feasible  solutions  is  a crucial  factor  in  the  performance  of 
the  branch  and  bound  algorithm.  Another  important  factor  is  the  branching  or 
cell  selection  rules  which  we  describe  next. 

(D)  Cell  Selection  (Branching)  Rules 

We  have  tested  several  of  many  possible  cell  selection  rules  and  have 
subsequently  reduced  our  choices  to  two  rules.  These  rules  utilize  both  the 
regret  concept  and  one  of  the  weak  lower  bounds,  WLBq,  discussed  earlier. 

In  order  to  describe  these  rules  we  first  define  the  following  quantities  at 
any  step  of  the  tree  search: 

F * the  set  of  columns  containing  more  than  one  basic  cell  having 
positive  flow 

U ■ the  set  of  columns  containing  no  fixed  in  cells 

B+  » the  set  of  basic  cells  having  positive  flow 

WLBq^ * the  weak  lower  bound  in  direction  x_  * 0 for  cell  (i,j). 

Using  this  notation  we  describe  two  branching  rules: 

(1)  BC- (basic  cell)  regret  rule 

a)  Column  selection  - choose  jeF  such  that  Reg(j)  > Reg  (k) 
for  all  kcF,  that  is,  choose  a column  having  largest  regret. 

b)  Row  selection  - find  the  two  smallest  cost  cells  (i^,j), 

(ij.J)  « B+  and  choose  cell  (i,j)  where  i - i^  if 

WLBQi  > WLBq^  , and  otherwise  i = i^;  that  is,  of  the 

two  smallest  cost  cells  in  column  j which  have  positive 


flows,  pick  the  cell  whose  weak  lower  bound  is  largest. 
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2)  G- (General)  regret  rule 

a)  Column  selection  - Choose  jeU  such  that  Reg(j)  > Reg(k) 
for  all  kcU.  (This  rule  is  the  same  as  (1)  (a).) 

b)  Row  selection  - find  the  two  smallest  cost  cells  (i^.j), 

(ijtj)  * I and  choose  (i,J)  where  i * i^  if 

WLBq^  > WLB^  , and  i ■ i^  otherwise;  that  is,  of 

the  two  smallest  cost  cells  in  column  j,  pick  the  cell 
whose  weak  lower  bound  is  largest. 

It  is  clear  that  by  using  the  G-regret  rule  there  are  more  cells  to  choose 
than  by  using  the  BC-regret  rule.  Also  the  BC-regret  rule  only  considers 
cells  (i,J)  whose  flow  x^  satisfies  0 < < b^;  we  will  call  this 

a fractional  flow,  in  agreement  with  the  usual  integer  programming  terminology. 
On  the  other  hand,  the  G-regret  rule  may  select  a cell  whose  flow  is  b ^ , i.e., 
an  integral  flow. 

In  most  integer  programming  algorithms  the  cell  selection  criteria  is 
limited  to  a choice  among  fractional  flow  variables.  The  reason  is  that  if  a 
non-constrained  variable  is  "naturally"  integer  then  it  would  seem  to  be  un- 
productive (i.e.,  would  increase  the  size  of  the  search  tree)  to  explicitly 
restrict  the  variable  to  be  integer  by  branching  on  it.  However,  it  is 
plausible  that  some  of  the  variables  which  are  naturally  integer  in  the  optimal 
solution  to  P might  also  have  the  same  integer  values  in  an  optimal  solution 
to  . Hence  it  might  be  worthwhile  to  restrict  these  variables  to  their 
respective  integer  values  in  hopes  that  the  size  of  the  search  tree  will  ulti- 
mately be  reduced.  In  line  with  this  reasoning  we  have  found  through  computa- 
tional experience  (section  6)  that  in  some  cases  the  size  of  the  search  tree 


i 


L 


is  smaller  when  using  Che  G-regret  rule  for  branching  rather  than  the 

BC-regret  rule  The  G-rule  usually  produces  a smaller  (in  number  of  nodes) 

search  tree  than  the  BC-rule  when  the  relaxation,  P,  has  a solution  which 

is  not  very  close  to  being  integer.  That  is  when  P contains  only  a few 

fractional  flow  cells,  the  BC-rule  concentrates  on  the  non-integral  (i.e.,  the 

non  single  source)  portion  of  the  relaxation,  and  in  doing  so  it  is  usually 

capable  of  performing  a quick  partial  enumeration.  However,  when  the  problems 

are  "tight,"  that  is  Z a - Z b is  small,  and  there  are  several 

i*I  jeJ  J 

fractional  variables,  it  is  often  better  to  fix  some  of  the  variables  (say  those 
variables  in  columns  having  a large  regret)  into  the  basis  even  if  they  are  not 
fractional.  This  avoids  the  possibility  of  these  variables  becoming  fractional 
as  we  move  deeper  into  the  search  tree.  This  is  accomplished  by  using  the  G- 
regret  rule.  Clearly  fixing  non-fractional  flow  variables  is  not  always  the 
best  alternative.  However,  we  have  noticed  dramatic  differences  in  the  per- 
formance of  the  algorithm  between  the  BC-rule  and  the  G-rule  on  some  of  the 
more  difficult  problems  which  we  have  generated.  For  example  on  a particular 
5 x 50  test  problem  (i.e.,  m ■ 5,  n ■ 50)  the  search  tree  using  the  BC-rule 
had  13,814  nodes  whereas  the  search  tree  using  the  G-rule  had  only  356  nodes. 

In  our  experience  the  BC-rule  usually  yields  a search  tree  with  fewer  nodes  than 
does  the  G-rule;  however,  the  BC-rule  is  not  invariably  better  than  the  G-rule. 
(E)  Statement  of  the  Algorithm 

Before  we  state  the  algorithm  we  define  some  notations: 
t ■ level  of  the  tree 

m the  transportation  problem  objective  function  value  at  level  L 
X.  » the  transporation  problem  solution  at  level  K 
U ■ the  set  of  columns  containing  no  cells  which  are  fixed  in 
T (K)  - (i,j)  if  cell  (i,j)  is  fixed  in  at  level  l,  and  * (i,j)+ 
cell  (i,j)  is  fixed  out  at  level  L 
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(Al)  The  Single  Source  Cost  Operator  Algorithm 

Step  (1)  (Heuristics)  Run  the  regret  and  the  largest  demand  heuristics  to 

get  an  upper  bound,  Z^.  If  neither  gives  a feasible  solution  then 

set  Z * 
u 

Step  (2)  (Initialize)  Let  U - J;  4-0.  Solve  P.  If  ZQ  > Z^  - 1 go 
to  (10).  Otherwise  go  to  (3). 

Step  (3)  (Variable  fixing)  For  all  non-basic  cells  (i,j),  if 

(c^  - (u^  + Vj))  * bj  > Z^  - ZQ  - 1 then  let  c^  ■ M otherwise 
go  to  (4). 

Step  (4)  (Cell  Selection)  Use  either  the  BC-regret  rule  or  the  G-regret 

rule  to  choose  a cell  (i,j)  upon  which  to  branch.  Save  WLBfl  . 

°i 

Step  (5)  (Fix  in  cell)  Replace  4 by  4+1.  Let  T(4)  » (i,j)*.  Fix  in  cell 
(i,J).  If  Z^  > (Z^  - 1)  go  to  (6).  Otherwise  if  X^  is  single 
source,  save  X^;  let  Z^  ■ Z^J  go  to  (6).  Otherwise  replace  U by 


U * 

■ { j ] and  go  to 

(4). 

Step 

(6) 

(Free  after  fixing  in) 

Free  cell  (i 

.J); 

replace 

4 by  4-1; 

go 

to  (7). 

Step 

(7) 

(Check  WLB)  If 

Z.  + WLBn  > Z go 
^ Wj  " u 

to 

(9).  Otherwise  go  to 

(8) 

Step 

(8) 

(Fix  out  cell) 

Replace 

J. 

4 by  4+1 . 

Let 

T (4)  - 

(i,j)+.  Fix 

out 

cell  (i,j).  If  > (Z^  - 1)  go  to  (9).  Otherwise  if  X^  is 
single  source,  save  X^J  let  Z^  = Z^;  go  to  (9).  Otherwise  go  to  (4). 
Step  (9)  (Free  after  fixing  out)  Replace  4 by  4-1.  If  4-0  then  go  to  (10). 
Otherwise  free  cell  (i,j).  If  T(4)  - (i,J)+  go  to  (9).  Otherwise 
go  to  (8). 

Step  (10)  Stop.  The  current  solution,  if  there  is  one,  is  optimal.  Otherwise 
there  is  no  feasible  solution  to  the  problem. 

(F ) Example 

Consider  the  problem  given  in  Figure  1.  Using  the  algorithm  just 
described  we  will  Indicate  the  various  steps  involved  in  finding  an  optimal 
single  source  solution  to  this  problem. 


Description  of  the  Operation 
Application  of  the  regret  heuristic  with  ■ 1 for 
all  jeJ  yields  Z^  ■ 391. 

U - {l,2,3, 4,5};  l - 0;  The  solution  to  P is  given 
in  figure  2(a)  Zq  ■ 279.7. 

Apply  the  non-basic  variable  fixing  rule  and  set 
c13  ” c24  * c43  “ C44  “ M- 

(using  Che  BC  regret  rule)  F » {2}  i^  « 1;  i^  » 2 . 
WLBq  - 1.41  x 16  - 22.6;  WLBQ  - 1.09  x 2 - 2.18. 


Choose  cell  (1,2). 

Let  c 12  - -M;  cn  - M;  c15  - M;  l - 1;  T(l)  - (1,2)'. 
is  a single  source  solution.  See  figure  2(b).  Save  X^. 


Let  c. 


.39,  cu  - 2,  c15  - 5.19,  l - 0. 


279.7  + 22.6  - 302.3  > 391. 
I - 1.  T(l)  - (1,2)+.  Let 


M.  ^ 


360.9.  X, 


is  not  a single  source  solution.  See  figure  2(c). 

F - {2,5}.  Reg (2)  - 111.6,  Reg(5)  - 37.9. 

Choose  column  2.  » 2,  “ 4.  WLBq  * 21.46  x 2 ® 42.9 

WLB-  ■ .04  x 16  * .64.  Choose  cell  (2,2). 

U4 

Let  c22  - -M,  c21  - M,  1-2,  T (2)  - (2,2)'. 

*s  single  source.  See  figure  2(d).  - 458. 


Let  c „ 


10.2,  C21  - 12.9,  c24  - 30.1;  0 - 1 


360.9  + 42.9  ■ 403.8  > 391. 


L m 0. 


Stop.  Z ■ 391  is  optimal,  x^  - 12,  x^2  ■ 18,  x,^  - 7, 

X34  " 10>  x45  ’ 26 • 


Figure  2(e)  contains  the  four  node  search  tree  for  the  example.  The 
first  number  in  the  parenthesis  at  each  node  is  the  objective  function  value 
of  the  current  relaxation  and  the  second  number  is  the  current  upper  bound  . At 
node  one  problem  P is  solved  and  the  solution  is  not  single  source  (see 
figure  2(a)).  The  upper  bound  at  this  point  is  371  which  was  obtained  by 
applying  the  regret  heuristic  with  cj  " 1 for  all  jej.  At  this  point  the 
non-basic  variable  fixing  rule  is  applied  to  fix  *24’  X43»  311,1  *44 

to  zero.  Cell  (1,2)  is  chosen  for  branching  by  the  BC  rule  and  at  node  2 
a single  source  solution  is  obtained  with  a value  equal  to  the  upper  bound. 

Next,  using  the  LIFO  rule,  we  move  to  node  3 where  the  solution  again  is  not 
single  source.  Then  we  choose  variable  (2,2)  and  move  to  node  4.  At  node  4 
the  value  of  the  relaxation  exceeds  the  current  upper  bound  thus  we  fathom 
the  tree  at  node  4.  Notice  that  the  weak  lower  bound  was  used  at  node  3 
so  that  it  was  not  necessary  to  fix  out  cell  (2,2). 

5.  SINGLE  SOURCE  VERSUS  ROW  UNIQUE  SOLUTIONS 

In  [ 7 ] the  concept  of  a row  unique  solution  was  used  to  characterize 
the  acceptable  solutions  to  the  sources  to  uses  problem.  In  the  present  paper 
we  have  replaced  the  idea  of  a row  unique  solution  by  a single  source  solution. 
Clearly  every  row  unique  solution  is  also  a single  source  solution,  but  the 
converse  is  not  true,  as  can  be  seen  in  the  example  shown  in  Figure  3(a).  The 
basic  solution  shown  there  is  a single  source  optimal  solution  (which  is  not  a 
row  unique  solution)  and  therefore  solves  the  single  source  transportation  problem. 
However,  the  reader  can  verify  that  the  primal  solution  shown  in  Figure  3(a)  is 
unique,  and  that  there  is  no  basic  optimal  row  unique  solution  to  the  problem. 

Hence,  to  find  a row  unique  optimal  solution,  it  is  necessary  to  generate  a 
branch  and  bound  search  tree  such  as  the  five  node  search  tree  shown  in  Figure  3(b). 
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This  tree  was  generated  by  using  the  cost  operator  algorithm  described  in 
section  4(E).  Because  the  search  tree  has  only  one  node  when  we  look  for  a 
single  source  solution,  and  five  nodes  when  we  look  for  a row  unique  solution, 
it  is  clear  that  use  of  the  single  source  solution  concept  yields  considerable 
computational  savings  even  on  this  small  problem  and  hence  much  larger  savings 
would  be  expected  on  larger  problems.  Thus  the  single  source  concept  is  es- 
sentially a way  of  circumventing  one  of  the  problems  created  by  primal  degeneracy. 


6 . COMPUTATIONAL  RESULTS 

In  this  section  we  will  discuss  the  computational  performance  of  the 
heuristics  (Sec.  3)  and  the  cost  operator  algorithm  (Sec.  4)  on  a set  of  randomly 
generated  problems  ranging  in  size  from  5 x 20  to  100  x 400.  These  problems 
were  generated  in  the  manner  described  in  [ 5]  as  follows.  We  use  a uniform 
probability  distribution  to  generate  random  integer  costs  c for  iel,  jeJ, 
between  1 and  50;  we  let  c^  “ 0 for  id;  similarly  we  generate 

random  integer  demands,  b^  for  j«J,  between  5 and  20.  Then  recalling  that 


is  the  smallest  entry  in  column  j,  we  set 


x.  . * b.  and  x.  * 0 for  i + j , all  jeJ 
jLJ  J ij  1 


We  then  calculate  the  largest  supply  , 


tsA ' .f.-u1 


which  is  needed  to  guarantee  feasibility  of  the  solution  in  (23). 


Finally  for  each  id,  let  the  supply  a,  ” QfS  where  the  slack  para- 


meter ot  is  chosen  to  be  less  than  1;  and  set  b 


I a - I b . 
id  jeJ  J 
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As  was  mentioned  previously,  > 0 is  a necessary  but  not  a suf- 

ficient condition  for  the  existence  of  a feasible  solution  to  P ' . In  general 
the  smaller  the  value  of  ot  and  hence  of  bn+j^»  in  a problem,  the  more 
difficult  it  is  to  solve.  When  the  slack,  , is  small,  there  is  more 

incentive  to  divide  the  flow  among  the  cells  in  a given  column,  and  thus 
violate  the  single  source  constraints.  Thus  by  varying  the  value  of  ot  we 
can  make  the  slack  for  a given  problem  larger  or  smaller.  For  at  * 1,  an 
optimal  solution  is  given  by  the  assignment  in  (23).  For  at  < 1,  the 
assignment  in  (23)  is  infeasible  and  thus  the  single  source  problem  is  likely 
to  be  nontrivial. 

Figure  4 contains  the  test  results  of  27  problems  ranging  in  size  from 
5 x 20  to  75  x 200.  All  of  these  problems  were  solved  with  a code  written 
in  Fortran  IV  using  the  BC-regret  rule  discussed  in  section  4(D)  on  a DEC-20 
time  sharing  system  at  CMU.  The  solution  times  in  Figures  4 and  5 are  subject 
to  some  measurement  error  due  to  variable  loads  on  the  time  sharing  system. 

A value  of  Cj  m (c^  * cjp/2  for  was  use<*  in  the  regret  heuristic. 

For  each  problem  the  regret  heuristic  chose  the  best  feasible  solution  out 
of  the  five  trial  solutions  which  it  generated.  We  decided  not  to  use  the 
largest  demand  heuristic  on  these  problems  since  the  regret  heuristic  always 
found  at  least  one  feasible  solution;  however  the  largest  demand  heuristic  may 
be  useful  in  other  problems  having  a different  data  structure. 

Notice  that  each  problem  size  is  repeated  three  times.  This  is  because 
the  difficulty  of  these  problems  can  vary  greatly  even  for  problems  having  the 
same  dimension.  Thus  we  felt  it  would  be  more  informative  to  solve  three  test 
problems  of  a given  size  in  order  to  demonstrate  the  potential  variation  in 
problem  difficulty. 
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In  addition  to  the  problems  shown  in  Figures  4 and  5 we  solved  three 
test  problems  given  to  us  by  Prof.  Terry  Ross.  Our  solution  times,  after 
taking  into  account  the  difference  in  computers,  are  comparable  if  not  slightly 
better  than  those  obtained  by  the  Ross  and  Soland  algorithm. 

The  regret  heuristic  error  was  evaluated  on  the  basis  of  a percent  error 
formula  which  is. 


percent  error  * 


X 100 


where  Zq^  is  the  optimal  value,  Z^  is  the  value  obtained  by  the  regret 


heuristic,  and 


Z„  = £ b.  min  {c  .}  (26) 

R jeJ  J is  I ij 

ZR  is  the  smallest  possible  objective  function  value  for  ?' . The 
purpose  in  subtracting  from  the  denominator  in  (25)  was  to  avoid  the 

problem  of  "scaling"  which  is  characteristic  of  network  problems  having  a 
transportation  structure.  That  is,  we  could  add  or  subtract  some  positive 
constant,  6,  to  each  of  the  costs  in  any  column  of  the  transportation  problem 
without  affecting  the  set  of  optimal  single  source  solutions.  Although  this 
scaling  does  not  affect  the  set  of  optimal  solutions,  it  does  however  affect 
the  value  of  an  optimal  solution  so  that  a standard  percent  error  formula 
such  as 


* Z0Pt' 


x ioo 


could  be  made  larger  or  smaller  by  scaling  the  data.  Thus  we  avoid  this  scaling 
problem  by  subtracting  Z_  from  the  denominator  in  (25).  For  example  if  6 is 

£v 


added  Co  the  cost  in  each  cell  in  column  j,  then  each  of  the  parameters  in 
(25)  increases  by  bj‘6  so  that  there  is  . ^ net  effect  on  the  percent  error. 
Also,  problems  for  which  Z q ^ = Z^  are  usually  uninteresting  and  thus  we 
attempt  to  avoid  this  situation  in  practice  by  making  O'  small  enough  so  that 
ZR  is  not  an  optimal  value.  A discussion  of  the  scaling  problem  for  general 
integer  programming  problems  is  given  in  Zemel  [10]. 

Note  that  the  CPU  time  to  run  the  heuristic  varied  approximately 
linearly  with  the  problem  size  as  m X n increases.  Note  also  that  the  percent 
error  of  the  heuristic  tends  to  decrease  as  m X n increases.  This  paradoxi- 
cal result  would  be  even  more  impressive  if  we  had  used  the  standard  error 


formula  (27)  instead  of  (25).  The  reason  that  the  heuristic  solution  method 


gets  better  as  m X n becomes  large  is  probably  due  to  the  fact  that  larger 
problems  tend  to  have  more  alternate  optimal  solutions,  and  there  is  the 
possibility  that  some  of  them  have  the  single  source  property.  In  any  case 
it  is  useful  to  know  that  a heuristically  generated  feasible  solution  to  a 
large  problem  has  a fairly  high  probability  of  being  optimal  in  those  cases  in 
which  it  is  not  possible  to  prove  optimality  in  a reasonable  length  of  time. 

The  variable  fixing  rule  (15)  seems  to  be  quite  effective  in  fixing 
non-basic  variables  to  zero.  In  many  of  the  problems  as  many  as  95  percent  of 
the  non-basic  bariables  are  fixed  at  the  initial  node  of  the  search  tree. 


As  we  mentioned  in  section  4(B),  if  the  variable  fixing  rule  consistently 
eliminates  more  than  say  607.  of  the  non-basic  variables  then  it  would  probably 
be  beneficial  to  use  a sparse  code  to  solve  the  transportation  relaxations. 

The  effectiveness  of  the  variable  fixing  rule  on  a given  problem  is  better 
the  closer  Z^  and  Z^  are  to  each  other,  because  then  rule  (15)  permits 
fixing  out  more  non-basic  variables. 
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We  tested  the  cost  operator  algorithm  both  with  and  without  the  weak 
lower  bound  calculations  and  found  that  the  weak  lower  bounds  are  indeed  cost 
effective.  They  help  by  reducing  the  size  of  the  search  tree  and  by  reducing 
the  overall  CPU  time.  As  can  be  seen  from  Figure  4,  the  weak  lower  bounds  can 
be  used  to  fathom  several  of  the  nodes  in  the  search  tree  which  would  other- 
wise have  required  explicit  enumeration.  Also  since  we  found  that  search 
trees  tend  to  be  smaller  when  the  weak  lower  bounds  are  used  for  branch 
selection  than  when  they  are  not,  we  conclude  that  they  provide  good  branch 
selections . 

As  expected,  the  single  source  transportation  problems  exhibit  a much 
higher  variance  in  their  execution  times  than  do  ordinary  transportation 
problems.  The  difficulty  of  a given  single  source  problem  depends  upon  many 
factors  such  as:  m,  n,  the  gap  between  Z^/  and  Z^,  the  amount  of  slack 
(i.e.  ’ the  cost,  supply,  and  demand  distributions,  the  problem  density, 

and  the  number  of  fractional  variables  in  the  initial  transportation  relaxation. 
Given  two  problems  of  the  same  size  the  one  having  the  larger  gap  is 
usually  more  difficult  to  solve.  Problems  in  which  the  demands  are  all  equal 
are  easy  to  solve  since,  as  we  mentioned  in  section  2,  partial  enumeration  is 
not  necessary  in  this  case.  If  all  of  the  supplies,  Si  iel,  are  greater 
than  or  equal  to  S,  as  defined  in  (24),  then  it  is  easy  to  see  that  the 
solution  given  by  (23)  is  optimal.  When  many  of  the  costs  are  very  large,  or 
equivalently  when  the  problem  is  sparse,  the  number  of  low  cost  solutions  is 
reduced  which  then  facilitates  the  implicit  enumeration.  We  have  found  that 
one  of  the  best  indicators  of  the  potential  difficulty  of  soving  a single 
source  transportation  problem,  ?' , is  the  number  of  fractional  variables  in 
the  transportation  relaxation  P.  A "large"  number  of  fractional  variables  in 
a basic  solution  to  P essentially  means  that  P is  not  a close  relaxation  to 
?'  and  thus  we  expect  that  such  problems  will  be  difficult  to  solve. 
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It  is  possible  to  determine  the  range  within  which  the  number  of  fractional 
variables  in  a basic  solution  to  P must  lie.  Consider  the  distribution 
of  the  basic  cells  in  any  ordinary  transportation  problem.  Given  an  m X n 
problem,  we  add  one  slack  column,  n+1,  to  the  problem  which  has  no  single 
source  restriction.  The  total  number  of  basic  cells  is  m + (n+1)  - 1 * m+n. 

Each  of  the  n+1  columns  must  contain  at  least  one  basic  cell  which  leaves 
m - 1 basic  cells  to  distribute  among  the  n+1  columns.  The  worst  case  (in 
terms  of  the  number  of  fractional  variables  in  the  basic  solution  to  P) 
occurs  when  m - 1 of  the  first  n columns  contain  two  positive  flow  basic 
cells  each,  in  which  case  there  are  2(m-l)  variables  which  violate  the 
single  source  criterion.  The  most  desirable  case  occurs  when  the  basic  solution 

to  P is  single  source,  that  is  when  all  of  the  m - 1 "extra"  basic  cells 

either  appear  in  column  n+1,  or  appear  among  the  first  n columns  and  have  a 

flow  of  zero  (they  are  primal  degenerate).  Thus  the  number  of  variables  which 
violate  the  single  source  criterion  lies  in  the  interval  (0,  2(m-l)).  Notice 
that  this  interval  is  independent  of  n.  Thus  we  can  increase  the  number  of 
column  locations  in  a given  single  source  problem  without  affecting  the  upper 
bound  on  the  number  of  fractional  variables.  In  practice  we  have  found  that 
when  a column  contains  more  than  one  basic  cell,  then  it  usually  contains  two, 
or  in  any  event  at  most  three  basic  cells.  Given  this  observation  one  would 


-27- 


it  a large  cost)  with  probability  .6.  Thus  approximately  60  percent  of  the 
cells  have  been  eliminated  from  consideration.  The  maximum  number  of  fractional 

cells  which  can  occur  in  any  of  these  problems  is  2 x 99  = 198.  The  number 

of  fractional  variables  encountered  ranged  from  13  to  46  and,  as  expected,  the 
problems  having  the  larger  number  of  fractional  variables  were  more  difficult 
to  solve.  This  is  also  true  in  figure  4 where  the  only  problem  which  could 

not  be  solved  within  10  minutes  of  CPU  times  was  a 75  x 200  problem  with  69 

fractional  variables  out  of  a possible  2 x 74  = 148.  Thus  we  believe  that 
number  of  fractional  variables  in  the  solution  to  P is  a good  measure  of  the 
potential  difficulty  of  a given  single  source  problem.  Utilizing  the  character- 
istics of  a class  of  problems  to  determine  how  difficult  they  may  be  to  solve 
is  important  since  as  we  mentioned  earlier  there  can  be  a large  variance  in 
the  difficulty  of  single  source  transportation  problems  of  the  same  size.  In 
figures  4 and  5 we  have  reported  the  characteristics  of  each  problem  which  we 
feel  provides  an  adequate  measure  of  their  complexity. 

If  we  look  at  the  total  CPU  times  in  figures  4 and  5 we  note  the  en- 
couraging fact  that,  except  for  the  one  75  x 200  problem  which  required  more 
than  600  seconds  to  complete,  the  overall  problem  difficulty  does  not  seem  to 
increase  with  increasing  problem  size.  (Of  course,  the  total  time,  which  in- 
cludes input  and  output  times,  does  increase  with  problem  size.)  Thus  it 
appears  that  the  single  source  transportation  model  can  be  useful  in  the  solution 
of  actual  problems,  especially  if  a user  is  willing  to  settle  for  a heuristic 
solution  whenever  total  CPU  running  time  becomes  excessive  . 


7.  CONCLUSIONS 

We  have  discussed  the  design  and  computational  testing  of  a cost  operator 
algorithm  for  solving  single  source  transportation  problems.  We  observed  that 
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the  computational  difficulty  of  a randomly  generated  test  problem  depended 
on  m,  n,  the  ?'  to  P gap,  the  amount  of  slack,  the  supplies,  the  demands, 
the  distributions  of  costs,  the  density  of  the  problem,  and  the  number  of 
fractional  variables  in  the  initial  solution  to  P. 

The  regret  heuristic  seemed  to  perform  very  well  on  these  test  problems. 
The  maximum  percent  error  was  41,  (as  calculated  in  formula  (25)),  and  more 
importantly  the  performance  of  the  regret  heuristic  improved  as  the  problem  size 
increased.  Using  the  value  of  the  regret  heuristic  as  an  initial  upper  bound, 
the  cost  operator  algorithm  solved  each  of  32  of  the  34  test  problems  in  less 
than  210  seconds  of  CPU  time  on  a DEC-20  computer. 

Two  surprising  observations  were  made:  First,  the  regret  heuristic  rule 
tends  to  make  smaller  errors  as  the  problem  size,  m X n,  increases.  Second, 
the  overall  solution  time  (with  one  or  two  exceptions)  does  not  tend  to  increase 
with  problem  size. 

These  two  observations  lead  us  to  conclude  that  the  single  source  trans- 
portation model  has  the  potential  of  becoming  an  easily  applied  operations 
research  tool,  even  though  it  is  used  to  solve  an  NP  complete  problem. 


1,  391)  optimal 
12“  18 


1) 

8 


Heuristics 


m X n 

1 

CPU* 

t 

i 

m 

e 

n 

0 

d 

e 

s 

CPU* 

c 

1 

m 

e 

Total 

P 

to 

p' 

Gap 

znrvr 

OPT 

ZR 

nodes 

fathomed 

by 

WLB 

Vars . 
fixed 
by 
rule 
(15) 

# of 

Fractional 

Variables 

Value 

of 

or  in 

Generator 

.05 

1.39 

85 

1881 

1730 

54 

46 

6 

.7 

Iml 

.04 

.16 

33 

2350 

2287 

2 

74 

2 

.7 

5X20 

0 

.04 

.16 

4 

1913 

1900 

1 

77 

2 

.7 

5X100 

36.9 

.49 

137.1 

28 

10847 

10533 

2133 

253 

6 

.7 

5X100 

9.9 

.48 

68.6 

28 

11732 

11421 

597 

340 

6 

.7 

5x100 

41.0 

.49 

19 

1.3 

5 

10487 

10409 

9 

365 

4 

.7 

5x200 

10.2 

1.56 

1766 

83.4 

14 

18769 

18573 

700 

728 

4 

.7 

5X200 

1.0 

1.60 

47 

4.54 

13 

19861 

19561 

15 

773 

6 

.7 

5X200 

20.7 

1.60 

2 

2.56 

0 

21538 

21258 

0 

692 

4 

.7 

40X40 

0 

.41 

7 

.97 

39 

762 

709 

3 

1425 

8 

.55 

40x40 

0 

.41 

23 

.98 

37 

1009 

913 

7 

1448 

4 

.55 

40X40 

0 

.46 

2013 

42.0 

81 

943 

818 

403 

1314 

20 

.55 

40X125 

11.5 

1.54 

177 

10.16 

9 

2523 

2497 

87 

4764 

24 

.55 

40X125 

0 

1.52 

2 

3.06 

6 

2770 

2757 

2 

4826 

4 

.55 

40X125 

0 

1.49 

173 

9.73 

5 

2365 

2360 

73 

4818 

16 

.55 

40X200 

14.2 

2.65 

27 

6.44 

2 

4150 

4115 

19 

7693 

16 

.55 

40X200 

0 

2.65 

76 

8.86 

9 

3803 

3784 

32 

7657 

17 

.55 

40X200 

3.5 

2.71 

3414 

209.2 

24 

4358 

4244 

1392 

7363 

28 

.55 

75X75 

5.1 

1.46 

861 

27.0 

29 

1172 

1133 

291 

5224 

20 

.50 

75X75 

7.5 

1.42 

804 

28.5 

38 

1269 

1203 

240 

5047 

33 

.50 

75X75 

0 

1.12 

0 

1.19 

0 

1158 

1158 

0 

0 

16 

.50 

75X125 

0 

2.16 

591 

28.8 

29 

1989 

1946 

231 

8738 

49 

.50 

75X125 

0 

2.02 

17 

4.65 

11 

2137 

2102 

7 

9051 

48 

.50 

75X125 

0 

2.00 

1 

3.77 

18 

1953 

1934 

1 

8916 

38 

.50 

75X200** 

mm 

4.34 

7153 

>600 

□ 

— 

3170 

1848 

13167 

69 

.45 

75X200 

SB 

4.12 

2 

8.5 

HI 

2958 

2942 

1 

14591 

33 

.45 

75X200 

H 

4.03 

15 

9.24 

El 

3009 

2979 

9 

14321 

35 

.45 

All  times  in  seconds  on  a DEC-20  computer. 

♦^oes  not  include  Input/Output. 

**  Terminated  after  10  minutes  CPU  time. 


Figure  4 


-32- 


m X n 

Heuristics 

n 

0 

d 

e 

s 

CPU 

t 

i 

m 

e 

Total 

P 

to 

p' 

GAP 

Z0PT 

ZR 

nodes 

fathomed 

by 

WLB 

Vars. 

fixed 

by 

rule 

(15) 

# of 

Fractional 

Variables 

Value 

of 

or  in 

Generator 

I 

CPU* 

t 

i 

m 

e 

100X100 

0 

2.11 

505 

23.88 

mOMMM 

2640 

159 

9221 

24 

.6 

100X150 

0 

3.30 

946 

48.07 

3251 

362 

14171 

30 

.50 

100X200 

3.3 

4.38 

772 

49.63 

4044 

3867 

290 

18543 

42 

.50 

100X250 

14.08 

5.73 

78 

16.02 

5987 

5916 

32 

24415 

13 

.55 

100X300 

0 

7.71 

20 

18.75 

6392 

6356 

8 

29386 

30 

.55 

100X350** 

- 

9.94 

6749 

>600 

— 

7036 

2751 

33486 

46 

.60 

100X400 

0 

10.78 

3 

23.39 

11 

8236 

8224 

3 

39292 

20 

.55 

All  times  in  seconds  on  a DEC-20  computer 
* > Does  not  include  Input/Output. 

**  Terminated  after  600  sec. 


Figure  5 


-33- 


RE  FERENCES 


[1]  Balachandran,  V. , "An  Integer  Generalized  Transportation  Model  for  Optimal 
Job  Assignment  in  Computer  Networks,"  Operations  Research.  Vol.  24,  No.  4, 
pp.  742-759,  July-August  1976. 

[2]  Christof ides , N.  and  S.  Eilon,  "The  Loading  Problem,"  Management  Sciences. 
17,  259-268,  1971. 

[3]  Demaio,  A.  and  C.  Roveda,  "An  All  Zero-One  Algorithm  for  a Class  of 
Transportation  Problems,"  Operations  Research.  Vol.  19  (1971),  pp. 1406-1418. 

[4]  Gaver,  D.  P.,  and  G.  L.  Thompson,  Programming  and  Probability  Models  in 
Operations  Research.  Brooks/Cole,  Monterey,  Calif.,  1973. 

[5]  Ross,  G.  T.  and  R.  M.  Soland,  "A  Branch  and  Bound  Algorithm  for  the 
Generalized  Assignment  Problem,"  Mathematical  Programming.  Vol.  8 (1975) 
pp.  91-105. 

[6]  Ross,  G.  T.  and  R.  M.  Soland,  "Modeling  Facility  Location  Problems  as 
Generalized  Assignment  Problems,"  Management  Science.  Vol.  24,  No.  3, 

Nov.  (1977)  pp.  345-357. 

[7]  Srinivasan,  V.  and  G.  L.  Thompson,  "An  Algorithm  for  Assigning  Uses  to 
Sources  in  a Special  Class  of  Transportation  Problems,"  Operations 
Research.  Vol.  21  (1973),  pp.  284-295. 

[8]  Srinivasan,  V.  and  G.  L.  Thompson,  "An  Operator  Theory  of  Parametric 
Programming  for  the  Transportation  Problem  - I,  II,"  Naval  Research 
Logistics  Quarterly,  19  (1972),  pp.  205-252. 

[9]  Srinivasan,  V.  and  G.  L.  Thompson,  "Solving  Scheduling  Problems  by 
Applying  Cost  Operators  to  Assignment  Models,"  in  Symposium  on  the  Theory 
of  Scheduling  (S.  Elmaghraby,  Ed.)  Berlin,  N.Y.,  Springer-Ver lag , (1973) 
pp.  399-425. 

[10]  Zemel,  Eitan,  "Functions  for  Measuring  the  Quality  of  Approximate  Solutions 
to  Zero-One  Programming  Problems,"  Discussion  Paper  No.  323,  Graduate 
School  of  Management,  Northwestern  University,  Evanston,  Illinois  60201, 

June  1978. 


Unclassified 

1 v tv. a • ivjm  OF  This 


mnwn  J««  citi«r««) 


REPORT  DOCUMENTATION  PAGE 

i''*~WC^OWT 


READ  INSTRUCTIONS 
□SFCRE  COMPl-ETTVC  "OHM 


I \ 
I 


y MSRR-429, 


. j.  covt  accession  moj  i.  *ecimickt,s  Cat  alog  numue* 

u t - </ ' /-  / ! 


/ j A SINGLE  SOURCE  TRANSPORTATION  ALGORITHM 


4 TITLE 


/ 


I i.  "»»e  of  mbmomt  4 »e;uoo  cc-.eReo  ■ 
Technical  Report  j 

February  1979  • 

! 1.  7CMF0MMING  OBG.  AePORT  nuMBCM 

I MSRR  429 

..  CONTMaC-'  ox  GRANT  NUMOifV.j 


I 


$1 


,?  lurnan'i/ 

fj  ( Robert  V./&agelhout 
Gerald  L./ Thompson 


(Js  j N00014-75-C-y)621  j 


1 


»exro«HiHa  organization  namc  amo  <oowu 
Graduate  School  of  Industrial  Administration 
Carnegie-Mellon  University 
Pittsburgh,  Pennsylvania  15213 


10.  AMOCHAM  liLEMEMT.  *«OJICT.  TAW 

AMCA  » WORK  OMIT  NUMBERS 

NR  047-048 


1 1 • 
I 


CSN-ac-UiMG  orrict  *««*  «no  acomcss 
Personnel  and  Training  Research  Programs 
Office  of  Naval  Research  (Code  434) 


S’ 


: nATa 

j j V Feb*BW|r*i79  / 

/pSfTBBWTM  0»  A>i£S 


32 


’ '4  uoNiTCMinG  -GENCY  name  A aOORCSKM  atlt—wml  trmm  c«amaIIIa«  OIUhi 

j u.  smcufiiTY  class,  a -**»»>«> 

j r 

\ iXii 

3Lp. 

! ! 

• j 

/ 

1 

l -ia.  050.  •S5I«ICATIu  ^uC-nMJau,NG 

! SCHEDULE 

1 ; 

r<7  it «i »u on  statement  ,«»  !*•*<>«) 

Approved  for  public  release;  distribution  unlimited. 


£ /,- 


y 


• ?.  OlSTAiauTtOM  STATCWlCMT  j/  cft«  •*•«■«<  w(«r«d  /«  IlMi  «0,  If  Sffl 


SUF#L*MCMTAAY  hotcs 


W KEY  *O«0S  (C«MHh 


i «i  titl#  l#  M4««Mry  *Y  * 


1 Transportation  Problem,  Single  Source  Transportation  Algorithm, 

Integer  Programming  Software,  Cost  Operator  Algorithm,  Generalized 
( Assignment  Problem. 


t 


,10.  ABSTRACT  'CaaiMm  aa  a.am  .!«•  I#  AAA  il—ultr  T A4»«A 

—A  A single  source  transportation  problem.  i.s.  an  ordinary  transportation  } 
problem  with  the  additional  requirement  that  the  entire  demand  at  each  demand  ; 
location  be  supplied  from  a single  supply  location.  It  is  a special  case  of  * 

! Ross  and  Soland's  generalized  assignment  problem.  Such  problems  occur  j 

i frequently  in  applications.  This  paper  gives  two  heuristic  solution  methods  ; 

. and  a branch  and  bound  algorithm  for  solving  single  source  transportation 
| problems.  A discussion  of  the  branching  rules,  variable  fixing  rules,  and  thej 
| computation  of  weak  lower  bounds  is 'given.  Computational  experience  wlth._t.hg_j 

V__  70;'m  _ ..  ...solution  of  randomly  generated  problems  having  up 

00  , JM  rj  |472  eoitiom  or  mov  «t  is  oriole . * t0  4Q.QQQ  integer  variables  is  reported. 

S/M  0 10  2"  0 I A*  660  I SECURITY  CL  AMI/ICATIOM  Of  THIS  RAOE  (Iha*  Of  /aiMSGj  ~ — 

Unclassified 


ft 


y 


/ > 


